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The  penetrator  problems  consisted  of  the  normal  impacts  of  rigid-body 
projectiles  into  sand  targets  at  a velocity  of  200  ft/sec.  The  projectiles 
for  Cases  1 and  2 were  6 in.  in  diameter  with  a 6-in. -long  ogive  nose 
(CRH  = 1.25).  The  Case  3 projectile  was  scaled  up  in  size  by  /3 , giving 
a 10.392-in.-  by  10. 392-in. -diameter  nose.  Two  penetrator  weights  were 
selected,  lAl.l  lb  (Case  l)  and  U2^.1  lb  (Cases  2 and  3),  to  give  projectile 
sectional  pressures  of  5 psi  for  Cases  1 and  3 and  15  psi  for  Case  2. 

The  solutions  of  each  penetrator  problem  were  run  from  initial  impact 
to  a penetration  depth  of  1-1/3  nose  lengths.  The  computed  penetration  dy- 
namics indicated  that  the  comparative  penetrator  decelerations  were  related 
to  the  sectional  pressure  by  a^tW/A)  - , in  agreement  with  conventional 
(rate-independent)  scaling  rules.  During  the  later  stages  of  the  penetrations, 
the  relative  decelerations  tended  to  depart  from  this  rule.  These  deviations 
are  believed  to  be  due  principally  to  the  effects  of  the  nonequivalent 
projectile  velocities  that  developed  during  the  penetrations,  resulting  from 
the  different  deceleration  levels  or  event  durations  associated  with  each  case. 

Since  there  is  some  experimental  evidence  suggesting  that  penetrator 
decelerations  of  this  type  are  related  by  aot(W/A)-l/2  , there  is  the  pos- 
sibility that  rate-dependent  behavior,  such  as  viscoplasticity,  may  be  occur- 
ring in  the  media  that  should  be  accounted  for  in  the  material  models. 

Detailed  results  of  the  numerical  analyses,  including  time  histories  of 
the  penetration  dynamics  variables,  force  loading  distributions,  and  field 
plots  of  the  target  response,  are  presented. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this 
verted  to  metric  (Si)  units  as  follows: 


Multiply 

feet  per  second 
inches 

kips  (force) 

pounds  (force)  per 
square  inch 


§1 

0.30U8 

0.025^ 

kkh8.222 

6.89k757 
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report  can  be  con- 

To  Obtain 

metres  per  second 

metres 

newtons 

kilopascals 


1 . SCOPE 


This  report  describes  the  results  of  a finite-difference  code 
calculation  of  projectile  penetration  conducted  by  California  Research 
and  Technology,  Inc.,  to  provide  information  of  use  in  the  study  of 
penetrator  scaling  relations.  This  calculation  was  for  the  third  case 
in  a series  of  three  penetration  analyses.  Pertinent  results  from  the 
analyses  of  the  first  two  cases’*-  are  included  herein  for  comparative 
purposes . 

2.  DESCRIPTION  OF  PENETRATOR  PROBLEMS 

The  penetrator  problem  considered  was  the  normal  impact  of  a rigid- 
body  projectile  into  a sand  target  at  200  ft/sec,*  as  shown  in  Figure  1. 
This  problem  was  designated  Case  3 in  a series  of  three  penetration  com- 
putations. The  conditions  for  all  three  cases  are  listed  in  Table  1. 
Cases  1 and  2 have  the  same  projectile  nose  geometry,  with  the  total 
projectile  weight,  or  W/A,  of  Case  2 being  three  times  greater  than  for 
Case  1.  The  Case  3 projectile  nose  geometry  is  scaled  up  by  a factor 
of  /3  over  that  of  the  Case  1 projectile,  but  retains  the  same  sec- 
tional pressure,  W/A  = 5 psi.  The  impact  velocity,  nose  shape,  and 
target  medium  were  the  same  in  all  three  cases. 

The  nominal  properties  of  the  sand  target  are  listed  in  Table  2. 

The  material  model  used  was  the  same  as  that  used  for  Layer  2 of  the 

Watching  Hill  Site  at  Defence  Research  Establishment  Suffield  in  a 
2 

previous  study,  except  that  the  low-pressure  end  of  the  yield  surface 
was  modified  so  that  the  cohesion  would  be  equal  to  zero. 

The  penetrator /target  interface  was  assumed  to  be  frictionless 
for  these  problems. 

* A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (Si)  units  is  found  on  page  3. 
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3.  COMPUTATIONAL  METHOD 


The  numerical  solutions  of  these  penetrator  problems  were  obtained 
with  WAVE-L,  a two-dimensional  finite-difference  Lagrangiar.  hydrodynamic- 
elastic-plastic  code.  This  code  has  been  employed  to  solve  penetrator 
problems  on  a number  of  Defense  Nuclear  Agency  and  U.  S.  Army  Engineer 
Waterways  Experiment  Station  projects  during  the  past  few  years. 

The  computational  grid  setup  for  the  Case  3 problem  is  partially 
shown  in  Figure  2.  The  basic  cell  size  used  in  the  target  was  ArQ  = 

Azq  = 0.866  in.,  or  1/6  the  projectile  radius.  Beyond  a radius  of 
10. U in.  and  a depth  of  20.8  in.,  the  cell  dimensions  were  gradually 
increased  in  2 percent  steps.  The  overall  grid  extended  to  a radius  of 
52.8  in.  and  a depth  of  63.2  in.,  comprising  2800  grid  points.  This 
grid  extends  beyond  the  boundaries  of  wave  propagation  that  took  place 
during  the  problem  (which  was  run  to  a penetration  depth  of  ^ lU  in. ) 
so  that  it  could  be  restarted  and  continued  to  longer  times  if  desired. 

The  integration  time  step  varies  during  the  course  of  the  problem, 
depending  on  the  sound  speed,  minimum  cell  dimension,  and  stability 
criteria  selected.  For  this  problem,  the  average  time  step  was  2.7  ysec. 

U.  NUMERICAL  SOLUTION  RESULTS 

Each  of  the  three  penetrator  problems  was  run  from  initial  impact 
to  a penetration  depth  of  1-1/3  nose  lengths,  i.e.,  penetration  depths 
of  8 in.  for  Cases  1 and  2 and  13*9  in.  for  Case  3.  The  computed  pene- 
trator deceleration  versus  depth  of  penetration  for  each  case  is  shown 
in  Figure  3.  Time  histories  of  the  penetrator  deceleration,  velocity, 
depth  of  penetration,  and  attached  length  (vertical  distance  from  nose 
tip  to  point  of  separation  of  projectile  from  target)  are  shown  in  Fig- 
ure U.  The  deceleration  for  Case  1 is  seen  to  be  approximately  three 
times  that  of  Case  2,  i.e.,  the  decelerations  are  approximately  in- 
versely proportional  to  the  penetrator  weights  (or  W/A).  Comparing 
Cases  3 and  1,  which  have  the  same  W/A,  the  peak  deceleration  levels  are 
seen  to  be  about  the  same,  with  the  deceleration  profiles  for  Case  3 
stretched  out  over  a greater  penetration  depth  and  time,  due  to  the  /3 
scale-up  in  size  of  the  projectile  nose.  Comparisons  of  the  computed 
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penet"ator  deceleration  and  the  change  in  projectile  velocity  during 
the  penetration  on  a scaled  basis  are  shown  in  Figures  5 and  6.  Using 
conventional  (rate-independent)  scaling  relations,  the  indicated  de- 
celeration and  velocity  change  ratios  should  be  equal  to  a value  of  1 
for  changes  in  just  the  weight  or  area  (diameter),  i.e.,  holding  all 
other  problem  variables  constant.  The  computed  results  are  seen  to 
vary  around  1 for  Case  3/Case  1.  For  Case  2/Case  1,  the  scaled  ratios 
are  seen  to  oscillate  around  a value  of  1 for  about  the  first  half  of 
nose  embedment ; the  ratios  then  steadily  increase  for  the  rest  of  the 
solution  (to  a projectile  embedment  of  1-1/3  nose  lengths).  The  gradu- 
ally larger  decelerating  force  seen  for  the  heavier  Case  2 projectile 
is  believed  to  be  due  to  a velocity  effect.  Note  in  Figure  3 that  the 
Case  2 penetration  velocity  is  higher  relative  to  Case  1,  with  the  gap 
increasing,  due  to  its  substantially  lower  deceleration  level.  The 

I 

relative  velocities,  plotted  versus  penetration  depth/nose  length,  are 
shown  in  Figure  7.  For  Case  3/Case  1,  the  deceleration  levels  are 
comparable,  but  the  event  time  is  scaled  up  for  the  larger  diameter 
Case  3 projectile,  so  that  the  relative  velocity,  vCase  3/V0ase  1’  :i's 
decreasing.  This  may  account  for  the  downward  trend  in  the  deceleration 
ratio,  aCase  3/acaoe  "beyond  one  nose  length. 

There  is  some  experimental  evidence  indicating  that  the  total 

1/2 

penetration  depth  is  proportional  to  (W/A)  for  low-velocity  projec- 
tile penetrations  into  soil.^’^  The  penetrator  decelerations  from  the 
jf  analyses  conducted  herein  do  not  deviate  largely  from  direct  inverse 

proportionality  on  W/A,  i.e.,  a<*(W/A)  1 , which  would  generally  in- 

dicate that  (for  these  type  penetrations)  the  total  penetration  depth 
would  be  directly  proportional  to  W/A,  i.e.,  Z^W/A  . The  direct  pro- 
portionality is  also  what  would  be  predicted  from  conventional  scaling. 
Conventional  scaling,  and  the  material  models  in  the  code,  presume  that 
rate  effects  are  not  significant.  The  difference  in  the  experimental 
and  computational  trends  suggests  the  possibility  that  rate  effects, 
such  as  viscoplastic  behavior  of  the  soil,  may  indeed  be  affecting  the 
dynamics  and  that  they  should  be  included  in  the  material  model. 

A series  of  field  plots  at  penetration  depths  of  2/3,  h/3,  6/3, 
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and  8/3  in.  were  obtained  from  the  Case  3 analysis,  as  shown  in  Fig- 
ures 8 through  15.  For  comparative  purposes,  these  depths  correspond  to 
the  same  penetration  depth/nose  length  values  for  which  plots  of  Cases  1 
and  2 were  obtained,  as  presented  in  Reference  1.  The  field  plots  show 
the  overall  projectile  and  target  configuration,  the  extent  of  plastic 
deformation  in  terms  of  the  generalized  plastic  strain  (e^),  the  veloc- 
ity vectors,  and  the  principal  stress  field. 

The  computed  loadings  on  the  projectile,  i.e.,  distributions  of 
normal  stress,  axial  force,  and  radial  force,  for  the  same  penetration 
depths,  are  shown  in  Figures  16  through  19. 
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TABLE  1.  PHNETRATOR  PROBLEM  CONDITIONS 


Projectile 

Target 

Impact 

Velocity 

Weight 

Diameter 

W/A 

Nose 

(tangent  ogive) 

Case  1 

200  ft/sec 

64.125  kg 

6 in. 

5 pel 

1.25  CRH 

Sand 

(141.4  lb) 

Case  2 

200  ft/ sec 

192.375  kg 

6 in. 

15  psi 

1.25  CRH 

Sand 

(424.1  lb) 

Case  3 

200  ft/sec 

192.375  kg 

10.392  in. 

5 psi 

1.25  CRH 

Sand 

(424.1  lb) 

(6  yT  in. ) 

1 

TABLE  2.  PROPERTIES  OF 

SAND  TARGET 

Density,  pQ 

1.42S6  gm/cm3 

89  lb/ft3 

Air  fraction,  Va 

401 

Cohesion 

0 

Unconfined  compressive  strength 

0 

Unconfined  tensile  strength 

0 

Mises  limit,  Y = JlJz 

max  max 

2906  psi 

Initial  value  of: 

Bulk  modulus,  KQ 

3292  psi 

Shear  modulus,  GQ 

2469  psi 

Constrained  modulus,  MQ 

6S85  psi 

Young's  modulus,  E0 

5926  psi 

Poisson's  ratio,  vQ 

0.2 

Dilatational  velocity,  cQ 

585  ft/sec 

10.392"  diam 


Figure  1. 


Penetrator  Problem  Conditions. 
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Figure  6.  Ratio  of  Scaled  Change  in  Velocity  versus 
Penetration  Depth/Nose  Length. 
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Figure  8.  Grid  Configuration  and  Generalized  Plastic  Strain  Field  at 
2/3-  in.  Penetration  Depth,  Case  3. 
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Figure  9*  Velocity  Field  and  Principal  Stress  Field  at 
2/3-in.  Penetration  Depth,  Case  3. 
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Figure  10.  Grid  Configuration  and  Generalized  Plastic  Strain 
Field  at  U/3-in.  Penetration  Depth,  Case  3. 
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Figure  11.  Velocity  Field  and  Principal  Stress  Field 
at  h/3-in.  Penetration  Depth,  Case  3. 
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Velocity  Field  and  Principal  Stress  Field  at 
6/3-in.  Penetration  Depth,  Case  3. 
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Figure  l6.  Penetrator  Loading  Distributions  at  2/3-in.  Penetration  Depth,  Case 
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Figure  IT-  Penetrator  Loading  Distributions  at  h/3-in.  Penetration  Depth,  Case 
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Figure  18.  Penetrator  Loading  Distributions  at  6/3-in.  Penetration  Depth,  Case 
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Figure  19.  Penetrator  Loading  Distributions  at  8/3-in.  Penetration  Depth,  Case 
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In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a facsimile  catalog 
caTd  in  Library  of  Congress  MARC  form?*  is  reproduced 
below. 
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